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We examine the dynamics of inflation driven by multiple, interacting scalar fields and derive a 
multi field version of the Hubble slow roll expansion. We show that the properties of this expansion 
naturally generalize those of the single field case. We present an analogous hierarchy of slow roll 
parameters, and derive the system of "flow" equations that describes their evolution, and show that 
when this system is truncated at finite order, it can always be solved exactly. Lastly, we express the 
scalar and tensor perturbation spectra in terms of the slow roll parameters. 
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^ ! I. INTRODUCTION 

Scalar field driven inflation [| solves the traditional problems (horizon, flatness, structure) of early universe cos- 
mology. Fortuitously, the need for an almost scale invariant perturbation spectrum strongly suggests that inflation is 
well described by the slow roll approximation [3, |3| : t he parameters that describe the departures from scale invariance 
^f-^ • are also those that must be small if the slow roll approximation is to hold [H . In practice, we normally consider only 
the two lowest order slow roll parameters. However, from a theoretical perspective we recognize this is a truncated 
f^) i expansion and we are implicitly assuming that there are no large high-order derivatives in the dynamical system 
V") ' describing inflation. 

For a single field, the general slow roll expansion was written down by Liddle, Parsons and Barrow 0. 1 The energy 
density and its proxy the Hubble parameter, H, is a function of the slowly rolling field, which serves as the independent 
variable of the dynamical system. The slow roll expansion is expressed as a hierarchy of differential equations, involving 
progressively higher order derivatives of H(cf>). In many practical applications, this hierarchy is truncated at some 
finite order, and the properties of the system have been carefully studied in recent years. In particular, Kinney [6| 
showed that this system possesses a set of attractors. Easther and Kinney |7j built on this observation, using the 
slow roll hierarchy to generate random models of inflation, and this technique was later employed in the analysis of 
the first year of WMAP data 0,0 and more recently in [HI lllj. However, the randomly generated models are not 
associated with a measure so one cannot quantify which region of parameter space is favored. Subsequently, Liddle 
J> , pointed out that the truncated slow roll hierarchy can always be solved analytically 01 ■ Leach and Liddle recently 
used this solution to implement and test different reconstruction algorithms, and their results illustrate how the lack 
of a measure makes it very difficult to draw quantitative conclusions about the parameter space of slow roll inflation 
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The developments above all refer to models of inflation driven by a single scalar field. However, multi field models 
are common in the literature and have a number of attractive features. Firstly, there is no theoretical guarantee 
that any inflationary phase can be described in terms of a single evolving field: the primarily motivation for doing so 
is simply an application of Occam's Razor. Moreover, many phenomenologically promising models of inflation (e.g. 
El El) 

require two or more interacting fields. Lastly, recent developments in string theory have motivated the study 
of multi field models in the context of the string landscape [ltl IT7L fiH. |n| . 

The purpose of this paper is to generalize the Hubble slow roll formalism to models with multiple fields. Section 
[n] recaps the slow roll regime for the single field case and Section IlIII gives the multi field version, which is the main 
result of this paper. We see that almost all features of the single field case carry over to the multi field case. In 
particular, the exact solution found by Liddle for the truncated single field slow roll hierarchy generalizes to the multi 
field case, and we write down an exact algebraic expression for the potential which corresponds to the truncated 
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slow roll hierarchy at arbitrary order. Finally we express the scalar and tensor perturbation spectra in terms of the 
generalized slow roll parameters. 



II. SINGLE FIELD FORMALISM 



A. Background Evolution 



In scalar field driven inflationary models, the inflaton field is governed by the following equation of motion, 

4> + 3H + = 0. 



(1) 



The spacetime background obeys the Friedman equations 



8ttG 



and 



H 2 = 



a 8nG 



a 3 



(2) 



(3) 



where H = ^ is the Hubble parameter. Equations (J2J and can be combined to produce the Hamilton- Jacobi 
equation [2(j 



H 2 W 



1-3^) 



8ttG 



v{4>), 



where the slow roll parameter e is 0, l2l| 



1 / I dH 



4ttG \H 



(4) 



(5) 



Note one freque ntly sees the potential slow roll expansion which can be obtained directly from the functional form of 
the potential |ll22| 



Vv 



16ttG ^ V(cj)) J 

1 V"{<t>) 
8ttG V{4>) ' 

We can understand the physical significance of e by rewriting © as 

= # 2 [1 - e] , 

and we immediately see that e < 1 corresponds to a > 0, the necessary condition for inflation. 



(6) 
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B. Hierarchy of Parameters 

The duration of inflation is most frequently parameterized by N, the number of e-foldings by which the scale factor 
grows during the inflationary phase, so a = e N . As a matter of definition, we obtain N by integrating H, 



V = J Halt. 



(9) 



As written above, N measures the number of e-folds before the end of inflation, and increases as t decreases. We can 
write N as a function of & via 



H 



N = / ^d<j) = 2VttG 



(10) 
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We now write a differential equation for e(N), 



where 



nVG 



1 (§H 



H 



H 



(11) 
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Equation I) 1211 is the first /Zoiu equation [|| and a depends on the second derivative of H . Defining additional parameters 
1 



l A = 



4nG 



" (d<f,Hy 



we create an infinite system of coupled first order differential equations: 



da 
dN 
d n X 

~d~N 



"X 



e(5CT + 12e) + 2( 2 A) 
1 



-(n- l)fi+ (n-2)e 



n + l 



X. 



(13) 

(14) 
(15) 



Looking at (| 1 31) we see that the derivative of m A depends on four (or fewer) quantities: m A, m+1 A, cr, and e. These 
flow equations can be truncated to a system of finite order if for some n, n X = 0, for all N. Once we have specified 
the values of the "A at some point N we can then calculated the full inflationary evolution as a function of N and 
then write down an expression for V(<f>). 



III. MULTIPLE FIELDS 



A. Background Evolution 



We now want to extend the Hubble slow roll expansion to an inflationary model driven by M scalar fields, fa. The 
Friedman equations (121/ill become 



H z = 



8ttG 



A I 



a _ 8ttG 
a 3 



j=0 
M 



The single equation of motion turns into M equations 



3^ + ^M=0, 



(16) 
(17) 

(18) 



where an <j> is now being used as a shorthand for all M fa. To produce multi field analogs of © and © we follow 
the method of [2l]]. Firstly, the time derivative of (|l(jf) is 



dH ^SnG^dfa fffa dV(fa) \ 
dt 3 ^ dt V dt 2 dfa J 



(19) 



2 Note that a = 2( 1 X) - 4e 
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which, when combined with the corresponding equation of motion 1|18[) . produces 

~dH 



dt 



dt 



(20) 



In general a time derivative can be changed into a derivative with respect to the fields by 

M , 



d d<j)j d 



dt ^ dt 

3=1 



and so 



2 Y^ 9H d<j>j 

2— F)rh: dt 



dt 



3 = 1 -'J - i=1 

Since the fields are independent, corresponding terms can be set equal, 

I dH 



4vrG d<j)j 



(21) 



(22) 



(23) 



This defines the unique relationship between the temporal derivative of a given field and the derivative of the Hubble 
parameter with respect to that field. 

We use this result to rewrite the equation of motion as [20I l2l| 



H z 



1 M 



12ttG 



dH 



8nG 



V 



(24) 



or 



H 2 



1 1 

1 ~ IP 12ttG ^ 



M /QH ,2 



8vrG, 



-V. 



We define a new set of parameters Cj, such that 



1 1 dH 



AttG H dfa 

In practice we will frequently look at the sum of squares, which can be written as a dot product: 



A I 



" 1 1 fdH x 2 



i=l 
r 2 



4ttG H 2 



(25) 



(26) 



(27) 



In principle, we should denote this quantity e but in practice we use e to emphasize the connection with the single 
field slow roll expansion, which is recovered by setting M = 1. Substituting (|2T|) into (|2*5|l we see 



1 ' 



BvrG 



(28) 



which is identical in form to the single field analog although 4> and e are now understood to have their multi- 
component forms. Similiarly, we can rewrite (|17H as 



-)=H 2 [l-e\ 
a 



(29) 



confirming that e < 1 



> for the multi field case. 



5 



B. The Multi Field Slow Roll Hierarchy 



Now that we have a well-defined parameter e, we can build our slow roll hierarchy. We start with 

dH dH 1 



dN 



1 



4ttG 





(0H\ 


dN 





1 



(30) 



or in terms of ^ = £ ^ gf, and «"' 



We define our first evolution equation: 



AttG 



^\ d 

^ dN 90" 

.7=1 J 



9 /Sff 



1 <9# 1 
H ~ e d0~H 



(31) 



dN 



M 

E 



6j Ajj 



(32) 



where Ajj a higher order parameter defined by 



Additional parameters are a generalization of (|33|l : 



i 



4ttG 7 i? 



i a 2 i? 



4ttG 



if' 



n 

i=2 



d m+1 H 



(33) 



(34) 



Note that the product in (|34|l begins at 2 and not 0; the first two indices are not included in this product. 

Using these higher order parameters, we can then calculate a system of coupled first order differential equations, 



d m \ r 



M 



M 



dN 



m+lA - 



+ m A c 



EE- 

\i=l k=2 



A; 



where we've used the relationship dcj>i/dN by looking at 

.(«) H 



N 



Hdt 



4>i{e) ei 



(35) 



(36) 



where the third equality is justified by the definition of e, (|26fl and (|23|) . 

We now explore the properties of this hierarchy. First, when M = 1 we have a single field and the lower indices 
become redundant; we drop them by 



i \ m \ 



(37) 



It is then trivial to check that we then recover the single field flow equations. 

These parameters exhibit several useful symmetries. Looking at l|34|) we see that the A are symmetric in the first two 
indices and separately symmetric in every other index. We can further reduce the number of independent parameters 
by recognizing that interchange of one the first two indices, a*, with one of the other indices, ay, is equivalent to 
multiplying by a factor of e ai /e a .. This is easily verified by a straightforward calculation. Consequently, we find that 
any two m A can be obtained in terms of one another (and the e 2 ) if they possess the same set of numerical subscripts 
as each other. At root, this corresponds to an implicit assumption that H is a well behaved function of the <pi and 
that all mixed derivatives commute. The final ambiguity that arises in the multi field case is that we are always free 
to rotate or translate the fields in any way, provided that the kinetic terms retain their canonical form. This will 
necessarily redefine the A's, although the overall physics remains unchanged. This issue does not arise in the single 
field case, since there is only one field and it cannot be transformed in a nontrivial way without changing the form of 
the kinetic term. 
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C. Alternative Parameters 

In developing a previous formalism for perturbations in multiple field slow roll inflation, Nibbclink and van Trent 
introduce the following parameters [23, |24| : 



r-M ~ (38) 



H 



n-l 



(39) 



where we use the ordinary time convention of [23 and the notation of |2J| . We can quickly show that the lowest order 
parameter used in their expansion is the same as our e: 

-= i = _f i™ = v^i^=f 2 = (am 

where our substitutions come from (|36fl . The similarity here allows both treatments to reduce the Friedman equation 
to its convenient form (|28ll . 

Our analyses diverge at this point. Nibbelink and van Trent choose to differentiate with respect to a temporal 
parameter (either t p3l | or a more general temporal measure r |24jp. whereas with reconstruction in mind we choose 
to define the parameters with respect to derivatives of the field. Since they both refer to the same physical system, 
the expressions obtained by Nibbelink and van Trent must be equivalent to ours. However, they differ in their relative 
tractability and in the clear connection between our formalism and the single field slow roll expansion. 

IV. EVOLUTION OF PARAMETERS AND RECONSTRUCTION 

The infinite hierarchy of parameters introduced in Section [ill Bl provides a complete description of the evolution of 
e = ejCj. We now follow the development of the single field analysis and truncate the hierarchy at finite order. As 
before, we see that if all the m X vanish for some to, then all the higher order parameters are also zero. To make this 
precise, assume there is some k such that for all to > k 

m+1 A Q0 ... Qm+ 4 =0 = 0. (41) 

The derivatives of these parameters vanish, thanks to (|35[) . What remains is a closed system of first order differential 
equations. 

We can now reconstruct the inflationary potential, analogously to the procedure outlined in ^yj. We begin by 
inverting 

/ q \ r 1 i 

(42) 

This requires knowledge of H(N), e(N), and <fi(N). We can trivially calculate e from the flow equations, while to 
track the evolution of the field one must look only at (|23|l 





'1 - 1; 


H (sttg) 






3 



' £i. (43) 



dN V 4ttG 

Using identities introduced previously, the differential equation for the Hubble parameter is 

i=l vl i=l 

where our final unknown is the initial value of this parameter. For this, we use the same normalization as in ; the 
amplitude of density fluctuations is on the order of 10 -5 so 

p 2n V e 
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From this it might appear that we have reconstructed the full potential. However, in order to put this to use we 
would need to measure a large number of m A with considerable precision, so one cannot use this to reconstruct the 
inflationary dynamics in practice. 

In and Q the single field flow equations were evolved numerically to recover the full inflationary dynamics 
after the slow-roll parameters were specified at a specific initial time. Subsequently, Liddle pointed out that the flow 
equations can be solved analytically to produce H{<j)) and V((f>). This treatment yields the inflationary potential in 
terms of the m A, as evaluated at some initial time. We now show that a similar truncation works in the multi field 
case. Assume that m+1 A 



a ...a m+ i 



for all choices of a,. Via lj!£f|> and this corresponds to 
gm+i H 



0. 



(46) 



The flow equations guarantee that the above equation holds for all values of (f>i . This means that a Taylor expansion 
of H around the origin of field space will contain only a finite number of terms, and provide an exact solution to the 
truncated flow equations. We use this knowledge to write H as 



H= £ A ijk ...<j>\ 

i,j,k... 



(47) 



In Ij47(l the sum runs from i + j + k + . .. = to i+j + k + ... = m so that the highest order terms in the expansion 
have total order m. 

Given an expression for H 7 we can find the coefficients ^4r//c... by taking appropriate partial derivatives. For instance, 
for any number of fields H contains a constant term ^4ooo... ; which can be found by evaluating H at 4>i — 



H 



0i=O 



An 



Ho- 



rn 



The value of Hq was, of course determined by (|45() . 

The higher order co-efficients may be best understood via a specific example. For argument's sake, consider the 
value of A121 in a three field model: 



2Ai 



d 4 H 



21 



where the prefactor arises from the multiplicity in the partial derivatives. As we know 



3 A 



122:! 



V4ttG 
= (J- 

or when we evaluate at the origin, using (|48[) . 



/ 1 \ 3 1 fdH\ fdH\ d 4 H 



H 3 \d4>2 J \d(f>3 J <90i<902<9</>2<9</>3 



H dcj)i d(j>2dc/)2 dcf)3 ' 



A 12 i = A 



000" 



(4ttG) 2 3 A 1223 



(49) 

(50) 
(51) 

(52) 



2 e 2 e 3 

From these coefficients, we can solve for the explicit form of e{(j>i) by taking appropriate partial derivatives of H(cj>i) 



1 



4ttG J H 



m 2^ \dcp k 



fc=i 



(53) 



The full analytic form of the potential is a complicated combination of these expansion coefficients (inverting 



v(4>) 



8ttG 
3 

8^G 



H 2 W 



1-3^) 



^2 A ijk...z<f>\<l>2<p3 



ii,j,k...z 



l 



12ttG 



M 

E 

fc=l 



,Q k (r,i,j,k...z A ijk...*<t>i< 



(54) 



(55) 
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Equation (|55|l provides an algorithm for finding an infinite number of exact inflationary solutions. This is simply the 
multi field generalization of a process that has long been applied to the single field case [2(| H^, |2(J |2]j - there is 
of course no guarantee that the resulting solutions will correspond to well-motivated potentials. Moreover, as noted 
previously we can transform the fields in any way that leaves the form of the kinetic terms unchanged. In this case, 
the explicit form of H(cj>) and V(<f>) will change (since these functions are not invariant under these transformations) 
although the physics is obviously not affected. Consequently, if one wanted to implement a multi field reconstruction 
algorithm, one would need to take this ambiguity into account. 



V. COSMOLOGICAL PERTURBATIONS 



We now turn our attention to the power spectrum of scalar perturbations psl l8(il | 



2 M 

E 

i=l 



dN 



(56) 



Consistent with |30| | we assume that N is the number of e-foldings from a reference point after complete reheating. A 
full discussion of isocurvature perturbations is beyond the scope of the current paper, and we do not include them here. 
For a clear treatment discussion of thegeneration of isocurvature perturbations in multifield inflation see Bassett, 
Tsujiakawa and Wands' recent review |29| . 

This problem has been solved for a very general set of scalar fields which couple via their kinetic terms as well as 
through their potential |3Ctl3l|. However, since we have assumed canonical kinetic terms from the outset, we can find 
considerable simplifications to the general expressions - in technical terms we are assuming that the metric on the 
scalar field space is flat, or that h a b = 5 a i, in equation (45) of |30|. The spectral index n-n is given by 



I = 



d\nP n 
d In k 



or, again 



by MM, 



n n -l 



t H 



G- 



d 2 v 



1 1 \ dN 



H 2 



id4> b SttG V 



c dN dN 
>l 3 d<t>i d<pj 



(57) 



(58) 



where we take advantage of summation convention. To evaluate this for our specific models, we use two important 
substitutions: 



( H 2 

= -2e 



H' 

a 



and 



d 2 V 



M 

= 'iH 2 [e a Eb + 1 X a b] — H 2 [ X ^bj 1 A a j + 2 ^abj] 

3=1 



(59) 
(60) 

(61) 



Equation becommes 

n n -l = -(2 + 4M)e- 



1 



M 

E 





u 


1 A a b 












e a £6 _ 





3 o,6=l 

In the single field limit, m = 1, we see that, to first order, 

n n - 1 = -4e + 2 A 
which agrees with the approximation used in 0- 



M i\ i\ _i_ 2 \ 

■y-^ ^bj A, , -r A,,/,; 

vi=i 



(62) 



(63) 
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The second observational parameter of interest is the ration of the amplitudes of the tensor to scalar perturbations, 
r = T/S. For this we use the combined results: the power spectrum for gravitational waves can be expressed as |3(j 



H ,2 



P g =[^) (64) 



and the spectral index for gravitational waves is given by [221 13 

H 



n g = 2^ = -26, (65) 

where the second equality translates the result into the slow roll parameters introduced in this paper. 
The spectra obey a generalized consistency condition [30l l33| , 

^ < Kl (66) 

which generalizes the consistency condition for inflation driven by a single scalar field. This ratio is proportional to 
the scalar:tensor ratio, r. 



VI. DISCUSSION 



In this paper, we have developed the Hubble slow roll formalism for multi field models of inflation. The properties 
of the multi field expansion are all easily understood as generalizations of their single field counterparts [4(. In 
particular, we show that there is a hierarchy of differential ("flow") equations that describe the evolution of the 
slow roll parameters. If only a finite number of terms in the hierarchy are non-zero at a given instant then the flow 
equations ensure that this truncation is preserved as the universe evolves. 

In principle, one could generate a large set of multi field models and attempt to constrain the slow roll parameters 
with observational data, as was done for the single field case via Monte Carlo reconstruction [7| • We briefly examined 
this possibility, but it appears that the larger number of parameters at each order in the multi field slow roll expansion 
makes it unlikely we can put meaningful constraints on the multi field parameter space for any realistic dataset. 
Moreover, while we have assumed that our fields have minimal kinetic terms, we can still rotate the fields into 
one another, preserving the form of the kinetic terms but modifying the functional form of the potential, and any 
reconstruction algorithm would have to take this ambiguity into account. 

While Monte Carlo reconstruction may not be practical in the multi field setting, Liddle, Parsons and Barrow's 
systematic treatment of single field inflation in terns of the Hubble slow roll expansion has found a large range of 
applications in thorough studies of the dynamical system that underlies these models. We consequently expect that 
this multi field generalization is a theoretical tool with a number of useful applications. We have shown there that 
there is a strong correspondence between the two systems and that many features of the single field case generalize 
naturally to their multi field counterpart. In particular, the truncated multi field flow equations can be solved exactly, 
yielding an expression for the potential analogous to that obtained from Liddle's exact solution to the single field 
hierarchy |12| . If one wished to do so, this solution to the flow equations could be used to generate any number of 
exact multi field inflationary solutions. Moreover, we can write down expressions for the perturbation spectrum using 
the multi field slow roll parameters, and in future work we plan to use this formalism to examine the perturbation 
spectrum produced when the inflaton does not follow a single, well-defined path in the multi dimensional inflationary 
potential. 
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